1 



STATIC RESPONSE FUNCTION FOR LONGITUDINAL AND 
TRANSVERSE EXCITATIONS IN SUPERFLUID HELIUM 



F. Dalfovo and S. Stringari 
Dipartimento di Fisica, Universitd di Trento, and INFN, 1-38050 Povo, Italy 



Abstract. The sum rule formalism is used to evaluate rigorous hounds for the 



(N 



o 

m 

(N 



C\l. 

^ ; density and current static response functions in superfiuid '^He at zero temperature. 

CN ! Both lower and upper bounds are considered. The bounds are expressed in terms of 

'j^ . ground state properties (density and current correlation functions) and of the inter- 

^:. 

atomic potential. The results for the density static response significantly improve 

Q I the Feynman approximation and turn out to be close to the experimental (neutron 



scattering) data. A quantitative prediction for the transverse current response is 
given. The role of one-phonon and multi-particle excitations in the longitudinal 
and transverse channels is discussed. 



PACS numbers: 67.40.-w, 67.40.Db 



I. INTRODUCTION 

It is well known that the linear response of a many body system to an external 
probe (density, current,. . . ) contains relevant information on the dynamic correla- 
tions among particles [1] . If the probe is static and coupled to the density of the 
system, one has to deal with the static density response function x(q'), related to 
the dynamic structure function S{q,uj) through the equation 

In the low-q' limit this quantity is fixed by the well known compressibility sum rule 

- limx(Q) 



where A^, M and c are the particle number, particle mass and sound velocity re- 
spectively. In the opposite case, g ^ oo, one finds the free-particle limit 

4NM 
- lim x(q) = —2-7 • 

In superfluid ^e the integral of Eq.(l) can be extracted, with rather good accuracy, 
from neutron scattering experiments [2]. An important characteristic of x{q) is a 
pronounced peak in the region of the roton wave vectors. It refiects the strong 
interaction between particles which tends to produce solid-like correlations in the 
system. This feature of the static response function plays an important role in the 
context of density functional theories for inhomogeneous Bose systems [3,4]. 

The calculation of x{q) represents a challenging theoretical problem. In this 
work we use the sum rule formalism to provide rigorous lower and upper bounds 



to the static response function at zero temperature. The basic idea was aheady 
introduced by Hall and Feenberg [5], who obtained bounds to xil) using the so 
called Feynman approximation. The difference between the lower and the upper 
bounds was, however, too large to make their approach of practical use. More 
recently an improved lower bound has been proposed [6] with the help of additional 
sum rules. In the first part of this work we will derive and evaluate explicitly new 
lower and upper bounds for the static density response function. The ground state 
properties which enter the relevant sum rules, needed to calculate the bounds, are 
taken from Monte Carlo calculations [7,8]. The resulting difference between the new 
lower and upper bounds turns out to be relatively small and, thus, they allow for a 
rather precise estimate of the true static response in very good agreement with the 
available experimental results at zero pressure. 

In the second part of this work we apply the same formalism to investigate 
the response to current excitations. We give a quantitative estimate of the static 
response in the transverse channel, through the use of a new lower bound, and we 
compare the results with the longitudinal case. 



II. DENSITY-DENSITY RESPONSE FUNCTION 

2.1 General formalism 

The linear response function characterizes the behavior of a quantum many 
body system subject to a small external perturbation. For static density excitations. 



at zero temperature, one can write the total Hamiltonian as 



H{X)=H + Xpl , (2) 



where 



^ = E^+E^(|r^-r.i) (3) 

is the unperturbed Hamiltonian, A is the strength of the perturbing field, and 

N 

pi = E^^''""^ (4) 

is the usual density operator. The static response function is defined as 

X(,)=lim<^ , (5) 

where |A) is the ground state of H{X). At zero temperature, standard perturbation 
theory gives xiQ) ^^ terms of the eigenstates and eigenvalues of H as follows: 

n 

With the usual definition of the dynamic structure function 

S{q,u) = J2\Hpi\0)\'S{o^-oOno) (7) 

n 

and its moments 

/•OO 

m,{q)= {hu;rS{q,u;)dw = J2\Hpl,\0)\H^^nor , (8) 



du = -2m-i{q) , (9) 



Eq.(6) reads 








X{q) = 


2 rS{q,u) 

Jo huj 
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The dynamic structure function S{q,uj) contains the detailed structure of the el- 
ementary excitations of the system. It is directly accessible to inelastic neutron 
scattering experiments on a relevant range of wave vectors. In particular, the in- 
tegral in Eq.(9) can be estimated rather precisely from the observed spectrum; the 
factor l/w makes the integral rapidly convergent at high w, reducing the effects 
of the complex structure of multiphonon excitations. An accurate measurement of 
x{q) is available, at present, only at zero pressure [2]. The theoretical determination 
of x{q) is a much harder problem; a direct approach, in fact, implies the calculation 
of non-uniform perturbed states with high enough accuracy to extract the linear 
limit; first results on this line are now becoming available [9]. In this work we choose 
an alternative approach based on the use of the moments nip with p > 0, which, 
differently from m_i, can be expressed in terms of known ground state properties 
with the help of sum rules. 

By using the completeness relationship ^^ \'n){n\ = 1 in Eq.(8), one finds 

moiq) = I duj S{q,u) = (p^Pq) (10a) 

rrniq) = jdu hu S{q,u) = ^([p^, [H,p^]]) (106) 

m2{q) = Jdu {huf S{q,u) = {[pi H][H, p^]) (10c) 

ms{q) = Jdu (huf S{q,u) = ^{[[pIH], [H, [H, p^]]]) . (lOd) 

The above equations relate properties of the excitation spectrum to mean values on 
the ground state. They are known as sum rules and have been extensively used in 
the theory of Bose liquids [1,5,6]. The moment toq coincides with the static form 
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factor 

NS{q)^moiq) , (11) 

where A^ is the number of particles. S{q) is related to the radial distribution function 
g{r) by means of 

S{q)-l = jdv{g{r)-l)e^^-' . (12) 

Several ab initio calculations are available for the pair correlation function g{r) 
(see for example Refs. 7 and 8). The corresponding S{q) well agrees with the 
experimental static form factor [10]. 

The energy weighted sum rule mi is the model independent f-sum rule 

which follows from the particle number conservation [1]. It can be easily derived 
from Eq.(lOb) taking into account that the interatomic potential commutes with 
the density operator. 

The 1712 sum rule can be expressed in terms of the current correlation function 
using the continuity equation 

[p^,H]=hq-J^ , (14) 

with the current density operator given by 

If q is taken along z, one finds 

m2(Q)=nV(4q'/.q) • (16) 
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Eq.(16) shows that m2{q) is proportional to the longitudinal component of the 
current correlation function in q-space. The TO2 moment can be also written in the 
following form [5]: 



^2^2\2 ^4 9 



<-^<'»l|^ -^-<'> 



'm2{q) = N 
where D{q) is the so called kinetic structure function 



(17) 



D{q) = /"rfrirfr2Cos(Q(zi -Z2))V?V^p(2)(ri,r2;r;,r'2) |ri=r;,r2=r^ • (18) 

While S{q) is fixed by the diagonal components of the two-body density matrix, the 
kinetic structure function D{q) requires the knowledge of the non-diagonal compo- 
nents. At P = we can use the Path Integral Monte Carlo calculations of Pollock 
and Ceperley [8,11] for the current density correlations in order to evaluate D{q). 
In Fig. 1 we plot the resulting curve. At q smaller than about 1 A~^ the accuracy is 
poor, due to the finite-size box of the PIMC calculations. At large q the numerical 
results are consistent with the asymptotic limit [5] 



2M 



\im D{q) = -^{KE) , (19) 



where (KE) is the mean kinetic energy per particle in the ground state. 

Finally we note that the cubic energy weighted moment ms can be rather easily 
evaluated carrying out the commutators in Eq.(lOd). One finds [12] 



msiq) = N 



^ W^' + 1^^™^ + fi^ Z"^" <7(r)(l-cos(q.r))(q.V)V(r) 



(20) 
where po and V{r) are the particle density and the interatomic potential respec- 
tively. 



2.2 The Feynman approximation 

So far we have shown how the sum rules too, toi, 1112 and ms can be determined 
from known properties of the ground state. Now we use them to fix rigorous bounds 
to rri-i and, consequently, to the static response function xil)- We notice that 
S{q, uj) is a positive function and the inequality 

rdu^^^il + ahufyO (21) 

Jo ri^ 

holds for any real a. Using the definition of the moments nip, one has 

m_i > -(2amo + ct^mi) . (22) 

One can vary the parameter a to make the r.h.s. of eq.(22) maximum. This yields 

a = — mo/mi and 

™-.>(^^^ . (23) 

mi riUpiq) 

where 

r^,(,) = H^ (24) 

is the energy of the phonon-roton excitation branch in the so called Feynman ap- 
proximation. Equation (23) provides a first rigorous bound to the static response 
function at T = 0. Using Eqs. (11) and (13) the same inequality can be written in 
the form [5] 

m_i(g)>m:^i(g) = -^-^ . (25) 

n q'^ 

The quantity m^i{q) corresponds to the Feynman approximation to m-i. 

In a similar way one can find an upper bound to x{q)- The crucial point is that 
in superfluid ^He there are no excitations with energy lower than the energy huJo of 
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the phonon-roton branch. Thus, it is possible to write 



which imphes [5] 



'%^rfc.< r^du , (26) 

n.^ Jo i^^o 



m_i < - — . (27) 

TiuJo 



Precise measurements of fiUo are available [2,13], so that the upper bound (27) can 
be accurately estimated at several pressures. One notices that the two bounds (25) 
and (27) would collapse in the exact to_i if the excitation spectrum were exhausted 
by a single collective phonon-roton mode. In that case also the Feynman energy 
fiuip would coincide with the true phonon-roton energy fiUo- 

In Fig. 2 the lower (25) and upper (27) bounds, calculated at zero pressure, 
are plotted as dashed lines, together with the experimental data. Experiments are 
consistent with the theoretical bounds which, however, turn out to be quite far 
each other. The sizable difference between the two bounds is a measure of the role 
of multiphonon excitations and is a signature of the inadequacy of the Feynman 
approximation . 



2.3 New bounds for rri-i 

Equations (21) and (26) can be generalized in a natural way through the proper 
inclusion of additional sum rules. Let's begin with the inequality 

r du^^^il + afiu + ^n'u^fyQ , (28) 

Jo nuj 



valid for any real a and /3. As before, we can write Eq.(28) as a lower bound for 
TO_i and vary both a and /3. After a straightforward calculation one gets [6] 



m-x{q) > 



m^iiQ) 



l-A{q)/e{q) 



(29) 



where 



<<l) 



and 



?7ii Too mo 



TO2 mi 

A(q) = 

nil mo 



m2 mi 
mi mo 



(30) 



(31) 



The ratio A/e takes important contributions from multiphonon excitations, through 
the moments m2 and ms. As a consequence we expect a significant improvement 
with respect to the Feynman approximation (25). 

As concerns the upper bound we generalize Eq.(26) in the following way [14]: 



"d..**"'-' 



huj 



Jo ri'^o 



(32) 



or, equivalently. 



m_i < h 27 



mi 

flUo 



mo +7 



2 / ^2 



hujr 



nil 



(33) 



Minimization with respect to 7 yields 



m-i{q) < 



rriQ 

hiVo 



1 



mo 
mi 



mi 
mo 



hUJr 



ni2 
rrii 



huJr 



(34) 



Note that since both mi /mo and m2/mi differ from huo, due to the important 
role of multiphonon excitations, Eq.(34) yields a significant lowering with respect 
to the Feynman upper bound (27). In Fig. 2 the new bounds (29) and (34) at zero 
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pressure are plotted as solid lines. At small q the microscopic ingredients used in the 
present sum rule analysis (density and current correlation functions of the ground 
state) suffer from a lack of accuracy due to the finite-size of the cell for Monte 
Carlo simulations. For this reason we have not shown the curves for the upper 
and lower bounds below q ~ 1A~^. As in the case of the Feynman approximation 
the experimental data fulfil the theoretical bounds. But now the two bounds are 
much closer each other over all the relevant range of g's. This means that Eqs. 
(29) and (34) provide an estimate of the static response function close to the exact 
value. We stress again that the evaluation of the two bounds (29) and (34) involves 
only ground state properties and, consequently, is much simpler than the explicit 
ah initio calculation of the static response function. 



III. THE CURRENT-CURRENT RESPONSE 

3.1 Longitudinal current excitations 

In this section we rewrite the formalism of Section 2.1 for the current response 
function. The current operator has been already defined in Eq.(15). As before, one 
adds a small perturbation to the Hamiltonian of the system. The perturbation is 
now a vector field proportional to the current density operator. The response of the 
system is given by the current response tensor 

■(0|J^q|ri)(n|Jtq|0) (0|jtq|n)(n|J^q|0)' 



Xm^(?,^) = X] 



(35) 



hu — huJnO -\- ifj huj + huJnO + ^V 

The transverse and longitudinal components of the response tensor can be studied 
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separately. Let's begin with the longitudinal one. We take q along z and define the 
longitudinal response function 

Then we define the quantity 

T^q^u) = J2\in\JlM'^(^ - ^no) (37) 

n 

and its moments 

/•OO 

m^{q)= {nujrT'^{q,uj)duJ = J2\HJU\0)fi^^nor ■ (38) 

The longitudinal static response function is the to ^ limit of Eq.(36). Using the 
definition of T^ one has 

-X^(,) = 2|rfc.I^ . (39) 

Indeed the determination of the static response function (39) is trivial. The key 
point is the continuity equation (14) which connects the matrix elements of the 
longitudinal current with the ones of the density operator: 



One easily obtains [1,15] 



WU^) = ^{n\p\\Q) . (40) 



L/ „x o^L /„^ o"^l(^) ^ 



-X^{q) = 2mtAQ)=2^ = - , (41) 

where mi is the density f-sum rule (13). In the same way one finds 



m2{q) 



m^il) = (J^U) = -rff (42) 

n q^ 
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and 

where 7712 and rris are the sum rules (16) and (20). The simphcity of results (41) 
reflects the fact that the response to a static longitudinal probe is, actually, a 
flctitious problem, related to gauge invariance properties [1]. 

3.2 Transverse current excitations 

The response to transverse probes plays a crucial role in the theory of super- 
fluidity. Actually the q -^ limit of the transverse response function deflnes the 
normal (non superfluid) density [1,15] of the system. This limit was extracted in 
Ref. 8 at several temperatures through a Path Integral Monte Carlo calculation of 
the transverse current correlation function and the use of the fluctuation-dissipation 
theorem. In this section we provide a flrst estimate of x^ at flnite q by calculating 
a rigorous lower bound at zero temperature. 

Let's rewrite Eqs. (36-39) by replacing the superscript L with T and the z- 
component of the current with an arbitrary component orthogonal to q. Similarly 
to the longitudinal case one has 

m^Q) = {JUJ.c,) (44) 

and 

m{{q) = \{[[JU,HlJ,^]) . (45) 

The calculation of these two sum rules follows exactly the procedure used for the 
1712 and 7713 sum rules in the case of density excitations (see Eqs. (16) and (20)). 
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In particular, the moment rriQ is the Fourier transform of the transverse current 
correlation function. The mf moment can be evaluated by carrying out explicitly 
the commutators in Eq.(45). One finds 



mi (q) = N 



■{KE) + ^ jdr g{r){l - cos{qz))VlV (r) 



3M2 
The transverse static response function is given by 



(46) 



-X^{q)=2ml,=2Jdu'^^^ . (47) 

The quantity T'^{q,uj) is positive, so that the inequality (22) holds even for the 
transverse moments and, thus, 

-X-M>2<4W^ , (48) 

This is a rigorous lower bound, valid at zero temperature. We have explicitly 
evaluated the r.h.s. of Eq.(48) at zero temperature, taking ttIq from the transverse 
current correlation function calculated in Ref. 8. The resulting lower bound is shown 
in Fig. 3. Since the PIMC calculations of Ref. 8 are carried out at low but finite 
temperature (the lowest value is T = 1.18 K), the use of the corresponding rriQ 
in Eq.(48) is meaningful only for q much greater than kT/c; the curve in Fig. 3 
corresponds to values of q well above this limit. 

At zero temperature the function y^ [q] should vanish at g = because the 
transverse current operator cannot excite phonons, which, in this limit, are the 
dominant excitations (the system is entirely superfiuid). At higher q multiphonon 
processes take place and x^(q') no longer vanishes. The first contribution to the 
static response is expected to be in q^ . In the opposite case q -^ oo one approaches 
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the free particle limit and all the ratios nip+i/mp tend to the same energy h q^ /2M. 
This implies 

n„. ,-(,). 2 ta MW)! . ^.«(^ , ,49) 



g— »oo 9^00 ra 



1 



"•^ fi q^ 



where we have used the asymptotic values 



2N 
\im m^ = —-{KE) (50) 

5^00 SM 



and 



The asymptotic behavior of x^(q) is shown in Fig. 3 as a dot-dashed line. The po- 
sition of the maximum of the solid curve provides a characteristic coherence length 
for superfluidity [1]. As expected, it is of the same order as the roton wave vector. 
The height at the maximum measures the strength of the interaction between par- 
ticles. In fact in a free Bose gas at zero temperature the function x^il) would be 
identically zero. 

Clearly Eq.(48) provides only a lower bound for —X^il)- A rough estimate 
of the difference between this lower bound and the exact value of the transverse 
static response can be made using the following arguments. First we note that only 
multiparticle excitations affect the transverse response function, since no transverse 
current is carried by the elementary excitations (phonons, rotons) of the system. If 
one assumes that the average energy and the spreading of multiparticle excitations 
are the same in the longitudinal as in the transverse channels, one can estimate 
the relative difference between m'^i and {m^)'^ / rn^ by analysing the multiparticle 
contribution to the corresponding longitudinal sum rules. Such a contribution can 
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be explicitly extracted from the experimental data on S{q, uj) [2]. With this proce- 
dure we conclude that the lower bound (48) should underestimate the exact value 
of — X^(q') by about 30% in the roton region. 



IV. CONCLUSIONS 

In this work we have investigated lower and upper bounds for the static response 
function of superfluid "^He at zero temperature. In the case of the density response 
x{q) the new bounds improve significantly the Feynman approximation, yielding 
estimates in agreement with the experimental data. In the case of the current 
response function we have given a first rigorous bound to the static response x^(q') 
and made a direct comparison with the longitudinal response, stressing the different 
role of one-phonon and multiparticle excitations. 
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FIGURE CAPTIONS 



Fig. 1 Kinetic structure function extracted from Eq.(17) with m2{q) and S{q) taken 
from Monte Carlo calculations [7,8]. In the limit q ^ oo one finds the asymp- 
totic value 0.8 A~^, as in Eq.(19). 

Fig. 2 Static response function for density excitations. Open circles: experimental 
values [2]; dashed lines: upper and lower bounds in Feynman approximation; 
solid lines: upper and lower bounds given in Section 2.3. 

Fig. 3 Static response function for transverse current excitations. Solid line: lower 
bound (48); dot-dashed line: g — > cxd asymptotic curve for the exact response 
function (see Eq.(49)). 
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